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Abstract. - We study the equilibrium and dynamical properties of a spherical version of 
the frustrated Blume-Emery-Griffiths model at mean field level for attractive particle-particle 
coupling [K > 0). Beyond a second order transition line from a paramagnetic to a (replica 
symmetric) spin glass phase, the density-temperature phase diagram is characterized by a 
tricritical point from which, interestingly, a first order transition line starts with coexistence of 
the two phases. In the Langevin dynamics the paramagnetic/spin glass discontinuous transition 
line is found to be dependent on the initial density; close to this line, on the paramagnetic side, 
the correlation-response plot displays interrupted aging. 



Various mean field spin glass models have been introduced and investigated in the last 
years; recently, a few models have been introduced mi^jElEl which combine features of spin 
glasses and lattice gas, and thus allow to describe important physical quantities, such as 
density. The frustrated Blume-Emery-Griffiths (BEG) model 3,4 , that we discuss here, is a 
general framework of such kind. Essentially it consists of a lattice gas [rii — 0, 1) where the 
particles interact through their internal degree of freedom [si = ±1) via quenched random 
couplings and also through a potential depending on a nonrandom coupling constant, K. A 
complete mean field equilibrium solution of the frustrated BEG has been recently found 0] 
by numerically solving the equations of the full replica symmetry breaking scheme through a 
suitable technique of integration. The solution has given evidence of a new interesting behavior 
in the spin glass panorama, namely a coexistence line between a paramagnetic (P) and a spin 
glass {SG) phase; the dynamics along this line has not yet been investigated. In this letter 
we study a spherical version of the mean field frustrated BEG model in the attractive case 
{K > 0) that can be analytically solved in its equilibrium and Langevin dynamics. 

The equilibrium phase diagram we get essentially reproduces that of the Ising version 
(sj = ±l,ni = 0,1) of the model the plane density-temperature shows a second order 
transition line from a P to a (replica symmetric) SG phase, ending in a tricritical point; here, 
a first order transition line starts where coexistence of the two phases exists. The tricritical 
point lies on the low density branch of the P spinodal {K < 1, Fig. la), or on the side of the 
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high density branch {K > 1, Fig. lb); in this last case both a, P/P and a P/SG first order line 
are obtained. In the Langevin dynamics from a random initial condition (rapid quench from 
high temperature) the equilibrium correlation functions display a diverging relaxation time 
at the second order line. Furthermore, the P/SG (and, if present, the P/P) discontinuous 
transition line is found to be dependent on the initial density d{Q) (Fig. 3); this critical line 
does not correspond to any singularity of the equilibrium correlation functions, but rather to 
a diverging plateau displayed by the density function d{t), when its long time limit jumps 
from an asymptotic (equilibrium) value to the other (Fig. 2). One can recognize the glassy 
behavior from the two-time quantities which exhibit an aging regime in the whole glassy phase; 
furthermore, still in the nonglassy phase, close enough to the P/SG discontinuous transition 
line, one can explicitly observe the pattern of interrupted aging from the correlation-response 
diagram (Fig. 4). 

Our spherical version of the frustrated mean field BEG model is defined as follows. Let us 
consider the Ising version of the model, described by the Hamiltonian 

H = - ^ JiySiTiiSjUj - ^ ^ niiij - fi'^ni (1) 

i<j i<j i 

fi being the chemical potential and Jy quenched Gaussian couplings with zero mean and 
variance \ Jfj\ — l/N. We introduce two new spin fields sn, S2i through the change of variables 
Si = Sii, rii = {siiS2i + 1) /2. Then, we assume Sii,S2i as continuous variables, subjected to 
spherical constraints (X^i — ^ f^'' = li2). Hereafter we shall use the variables Sai to 
describe the model. The statics can be solved using the eigenvalue density p (A) of a large 
random matrix 5 . The P and the SG phase are described by: 



W " z + 2/3m + WKd (2) 

z = W d = 1 - — — — — — 9 = d - 4 SG 

^ 2/3 + 2/3^1 + 2f3Kd ^ /? 

where d = [(siiS2i) + 1] /2 denotes the density, q — {sn + S2i)^ /4 the overlap and z = 
z\ — — fiKd with z\ = zi the two Lagrange multipliers related to the spherical constraints. 
The stability analysis yields two critical lines, called spinodal lines, at which the corresponding 
solution becomes unstable against density fluctuations (and both the compressibility and the 
specific heat are found to diverge): 



-1 
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T= |^-ii:+yA'2 + (l-d) P-spinodal; T = 2i\: (1 - d)^ 5G-spinodal (3) 

A second order transition line is located by d = T in the plane d—T (or by T = Tc (/x, K) in the 
plane fi — T); below this line one finds from the first of Eqs. Q that the P solution becomes 
unphysical. By increasing the parameter K the spinodal lines move to higher temperatures, 
while the second order line is ii'-independent. For K < 1 (Fig. la) the second order line meets 
the spinodal lines in a tricritical point A, with temperature Ta, located on the low density 
branch of the P spinodal: 



2K-3-2^^+^i2K-3-2^lf + ^6K{^^ + l) 4K + 1 - VSKTT 
Tci^^,K) = — ; Ta{K) = — 

(4) 
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Below the point A, the system undergoes a P/ SG first order transition where the order param- 
eters q and d discontinuously jump and a latent heat is involved; the critical line is numerically 
determined by imposing the free energy balance between the two solutions. Between the spin- 
odal lines in the plane d — T, no pure phase is achievable, not even as a metastable one, and 
the system is a mixture of the two phases. In the zero K limit, Ta vanishes and the first order 
line disappears 7,. For K > 1 (Fig. lb) the second order line intersects the P spinodal in 
its high density branch; in this case a first order transition from the low to the high density 
P phase is also exhibited. Analytical expressions of the thcrmodynamical quantities can be 
derived in terms of the above solutions, Eqs. The linear susceptibility, the compressibility 
and the specific heat show a cusp by crossing the second order line 7 and a finite jump at 
the first order ones. The SG susceptibility is infinite in the whole SG phase, while remains 
finite in the P phase, diverging when the second order line is approached [7]- 
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Fig. 1 — Equilibrium phase diagrams in the plane fJ. — T and d — T for K = 0.5 (a) and X = 3 (6). One can 
have a spin glass (SG) phase or a paramagnetic (P) phase (which forks at the point iJ in a low density P~ 
and an high density P+ component). The transition lines are marked in bold; the spinodal lines, at which 
the relative phase becomes unstable, are also shown as continuous lines; dotted lines represent unphysical 
analytical continuations of critical lines. Finally, full circles mark the tricritical point A and the bifurcation 
point B on the P spinodal line. Notice that the P phase becomes unphysical below the second order line in 
the plane d — T; a P/ SG first order transition line starts out from the tricritical point A, located on the P~ 
spinodal in the case (a) or on the side of the P"*" spinodal in the case (b). 



Now we deal with the dynamics. We assume that the two spin fields Sai evolve via usual 
Langevin equations; using the mean field approximation we can replace the particle-particle 
interaction with —Kd{t) {siiS2i + 1) /2. Thus when projected onto the basis A where Jy 
is diagonal, the Langevin equations become: 

- r^i^ ^ (^/^ " ^1 /2 A/4 +{n + Kd (t)) /2\ /sia\ AiA (t) + CiA {t)\ 

dt \s2x) Va/4 +{^l + Kd [t)) 12 A/4 - Z2 {t) /2 ) \s2x) ^ \h2x (t) + 6a (t)) ^ ' 

where Za{t) (a — 1,2) are two Lagrange multipliers enforcing the spherical constraints, 
haX (t) two external fields and (t) the thermal noises with zero mean and correlations 
(CaA (i) Cfe/i = '^T6ab5\^5 [t — t'). The initial fields s^a (0) are assumed to be Gaussian 
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variables with zero mean and variance Sa\ (0) Sb\ (0) = 1 + 2 (1 — Sat) {d (0) — 1) where d (0) is 
the initial density. We study the correlation functions Cia {t^t') — J dXp{X) {si\ (t) Sa\ {t')) 
(a = 1, 2), the response functions Gia {t, t') = J dXp (A) 6 {si\ (t)) /5ha\ (f) and the density 
d (t) = (Ci2 {t, t) + 1) /2. Due to its linearity, the Langevin system Q can be explicitly solved 
for given noises, thus the dynamical quantities can be evaluated averaging over the noises and 
the eigenvalues A 7 , 8 . We find for t > t' and zero external fields: 



Cla {t,t') = 



d\'-±^\-2T 



+'na 



1 



Gla {t,t') 



't + t' 



I t it') hit ~t') l ^jt') 
v{t) t-t' 2 y E (t) 



(*+*')/2 h{t + t'-2u) V{u) 

t, t + t'-2u r(^±^) 

;(*+*')/2 



(6) 



where -qa = {-if^^-, h [t) is the modified Bessel function of order 1 and S (t) = T [t) e^^^* d(")d«+2Mt^ 
while r [t) and d{t) satisfy the following coupled integro-differential equations: 



T{t)d{t) =d(0) 
d 



h {2t) 



T 



dt 



d{t) 



1 







t 



u 



dlt)) —liir{t) + 2fi + 2Kd{t) 
dt 



(7) 



T 



with the initial conditions F (0) = 1 and arbitrary d (0). 

The special case if = is analytically tractable by standard Laplace techniques [7] and one 
can find the full time dependence of F and d, but this is not the case for K > Q. However, it 
can be verified that the long time fixed points of the d (i)-equation are described by the P and 
SG equilibrium equations (O with F (i) respectively proportional to e*^^/'^^* and e'^'^ / V^Trt^ ; 
furthermore local analysis of these fixed points yields the same stability conditions of the 
equilibrium theory. Moving from high temperature towards the equilibrium second order line, 
a diverging equilibration time t^^'^ {T — Tc)^"^ is foimd 7 . In the (both P/SG and P/P) 
first order region, between the spinodal lines in the plane /i — T, two distinct asymptotic 
solutions are possible; to know which of them is established for long times, stability analysis 
only is not sufficient, but an appropriate matching with the initial conditions is required. This 
has been achieved by numerically solving Eqs. Q). As an example, the different curves in 
Fig. 2 are obtained by keeping d (0) , K, T fixed and taking several values of /i slightly above 
and below the estimated critical threshold fjid at which the system jumps from the P to the 
SG solution. The approach to the critical threshold of the density function is accompanied 
by a diverging plateau [rlf ~ — In {p,d — m)) at an unstable fixed point given by the further 
root at the SG of Eqs. before falling in one of the two stable asymptotic values. Fig. 
3 shows that at fixed K the whole critical line depends on the value of the initial density 
d(0). Lower values of c?(0) favor the P phase and, for very low values, the SG solution is 
obtained only when the P solution becomes unstable, i.e. at the P spinodal line; in this last 
case Tg^* ~ (fip-sp — fJ-)^^- A such behavior is to be essentially related to the competition of 
the two phases in the first order region; the same features are found for the P/P discontinuous 
transition, apart from the power law behavior of F (t) e"^* visible in Fig. 2. 

The equilibrium dynamics in the P phase is formally equivalent to a trivial schematic 
formulation of the Mode Coupling Theory (MCT) with two modes a = I, 2 (using the notation 
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of |H1 the MCT kernel is Fi); the two-time quantities 10 are functions of t — t' and satisfy 
the fluctuation-dissipation theorem (FDT), decaying to zero with relaxation time rf^"^ [7]. 
In the SG phase, due to the power law factor in F (i), a nonequilibrium long-time dynamics 
is found; this is characterized by a regime {t ~ t') where FDT holds, followed by an aging 
regime (t > t') where the correlation Cia is a function of t' /t and the integrated response 
Xia (^j = /(' Gia {t, u) du plotted as a function of the correlation, xia (Cia), is constant |7I8| . 
We notice that a SG-\ike nonequilibrium regime is present also in the P phase, close enough 
to the P/SG discontinuous transition line, for waiting times lower than r^^* (Fig. 2). In this 
case one can well distinguish the following three steps in the relaxation process (see the curves 
Txii (Cii) in Fig- 4): an initial FDT regime (t^^* > t ~ t') followed by an aging regime 
i'^eq^ > < > t'), as for the SG phase; finally an interrupted aging regime (t > r^^* > t'). 

In conclusion, we have analyzed the equilibrium and dynamical properties of the spherical 
frustrated BEG model at mean field level for > 0. The introduction of dilution and 
attraction in a standard spin glass model generates a coexistence line between the P and 
SG phase that exhibits new interesting features. Future extensions or improvements of our 
framework can be realized, such as a version of the model with p > 2-particles interactions |10| ; 
this last case could display interesting connections with recent studies on the MCT for glassy 
systems with attractive interactions 
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Fig. 2 — Graphs of density d it) and F (t) e~^' (in logarithmic scale) when the system is close to the discon- 
tinuous P/SG transition. The curves are obtained by numerically solving Eqs. for M ~ A'd = ±10~^ with 
a; = 2, 4, 6, 8 from left to right and ~ —1.27; the other parameters are fixed sA, K = 1,T = 0.2 and d (0) = 1. 
Before falling in one of the two stable asymptotic values the density function exhibits a plateau at an unstable 
intermediate fixed point, whose length diverges as approaches the critical threshold /i^j (t^^* ~ — In (/i^j — p.)). 

Fig. 3 — We show the P/SG discontinuous critical lines for three values of the initial density, d(0) = 1, 0.75, 
0.5 (squares, circles, triangles), in the region of the plane fJ. — T between the spinodal lines [K = 1). The 
equilibrium first order line is also shown (dashed line). For d(0) lower than ~ 0.4 the critical line coincides 
with the P spinodal line. 

Fig. 4 — We plot Txii vs. Gn for several values of t' (indicated in the figure) when the system is in the 
P phase but very close to the P/SG discontinuous transition line (K = 1, T = 0.1, d{0) = l,fj, ~ —1.44) 
and t' < T^f ~ 100. One can distinguish in the relaxation process an initial FDT regime (Txii = 1 — Cn), 
followed by an aging regime, as for the SG phase, where the curves approach the dotted line {Txil = 1 + T — d 
where d is the density at the unstable fixed point); finally an interrupted aging regime, due to the finite 
equilibration time, in which the curves bend upward. 
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